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Abstract 

The  Naval  Research  Laboratory  (NRL),  as  part  of  the  Global  Positioning  System  (GPS) 
Center  of  Expertise  (COE)  test  program,  performs  absolute  calibration  of  GPS  receivers  for 
time  transfer.  NRL  did  a  simulation  of  this  calibration  method  to  verify  the  procedure.  The 
simulation  focused  on  the  effects  of  filters,  external  to  receiver,  on  receiver  calibration.  Two 
areas  were  investigated  in  this  simulation.  Will  changes  in  the  data  from  two  different  types  of 
receivers  track  if  the  bandwidth  of  the  signal  source  to  both  receivers  change?  If  a  filter  is 
added  between  the  antenna  and  a  receiver,  can  this  combination  be  calibrated  in  parts  or  must 
this  combination  be  calibrated  together? 


CALIBRATION  METHOD 

The  calibration  of  a  GPS  receiver  is  accomplished  with  a  GPS  signal  simulator.  The  simulator  generates 
two  signals:  a  simulated  GPS  antenna  signal  and  a  one  pulse  per  second  reference  [1],  The  receiver  is 
connected  to  the  simulated  antenna  signal  and  the  time  determined  by  the  receiver  is  compared  with  the 
time  of  the  simulator  reference  (Figure  1).  An  additional  step  of  calibrating  the  simulator  is  also 
necessary,  which  is  done  by  connecting  a  high-speed  digital  scope  to  both  outputs  of  the  signal  simulator. 
The  pattern  on  the  scope  is  shown  in  Figure  2.  The  calibration  of  the  simulator  is  defined  as  the  delay  of 
the  null  in  the  envelope  of  the  GPS  signal  for  the  one  pulse  per  second  reference.  In  some  cases,  an 
antenna  system  is  placed  between  the  simulator  and  receiver.  The  combination  can  be  calibrated  as  a 
whole  or  in  parts.  If  in  parts,  the  delay  of  the  antenna  system  is  measured  with  a  network  analyzer.  The 
absolute  calibration  of  the  receiver  is  defined  as  the  delay  of  the  receiver  tick  from  the  simulator  tick  less 
the  simulator  calibration  and  less  the  antenna  system  delay  if  the  calibration  is  done  in  parts. 
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BACKGROUND 

A  simple  example  will  illustrate  the  concern  NRL  has  about  the  calibration  method.  The  example 
contains  a  single  pole  filter  placed  between  a  bi-phase  modulator  and  a  Delayed  Locked  Loop  (DLL) 
receiver.  The  output  of  the  single  pole  filter  will  have  an  exponential  envelope  after  a  bi-phase  transition. 
The  envelope  will  be  of  the  form 


Envelope  =  1 ,  if  t<0 
Envelope  =  (2e't/Tc  -1),  if  t>0 

where  t  =  0  at  the  time  of  the  transition.  The  time  constant,  Tc,  is  equal  to  twice  the  reciprocal  of  the 
bandwidth  of  the  filter  in  radians  per  second.  The  delay  of  the  filter,  as  measured  by  the  network 
analyzer,  will  be  equal  to  the  time  constant.  This  transition,  when  viewed  with  the  digital  scope,  will  have 
an  envelope  null  at  t  =  .693 1  Tc . 

The  DLL  will  try  to  track  this  transition.  A  DLL  has  an  early  and  a  late  correlator.  The  code  reference  to 
late  correlator  is  delayed  from  the  reference  to  the  early  correlator  by  a  time  called  the  correlator  spacing. 
This  is  less  than  or  equal  to  a  code  chip.  The  loop  error  signal  is  the  difference  in  the  output  of  the  two 
correlators  [2].  If  the  code  reference  to  the  correlators  is  not  changing,  the  difference  in  the  correlator 
outputs  will  be  zero.  The  only  time  the  loop  error  signal  is  nonzero  is  after  a  reference  change  for  a 
period  of  time  equal  to  the  correlator  spacing.  The  point  where  the  DLL  will  track  this  transition  can  be 
found  by  finding  where  the  cross-correlation  of  the  envelope  and  the  correlator  spacing  is  zero.  If  this 
spacing  is  very  small,  the  DLL  will  track  the  envelope  null.  Since  the  slope  of  the  envelope  is  greater 
before  the  null  than  after  the  null,  larger  correlator  spacing  will  cause  the  DLL  to  track  after  the  null.  The 
difference  between  the  envelope  null  and  the  DLL  tracking  point  is  a  function  of  the  filter  and  the 
correlator  spacing.  With  a  single  pole  filter  of  bandwidth  equal  to  twice  the  C/A  code  rate,  the  difference 
between  the  zero-crossing  point  and  the  one  time  constant  point  can  be  47.7  ns. 

The  three  measurements:  the  digital  scope,  the  network  analyzer  and  the  DLL  all  give  different  results. 
Figure  3  shows  a  decay  envelope,  the  one  time  constant  point,  and  where  the  DLL  would  track  if  its 
correlator  had  a  spacing  of  two  time  constants. 


SIMULATION 

The  methods  used  in  the  above  example  will  not  support  complex  filters.  A  Fast  Fourier  Transform 
(FFT)  method  was  selected  for  a  simulation  because  it  supports  complex  filters,  transforms  between  time 
and  frequency,  and  correlation.  Figure  4  is  a  block  diagram  of  the  processing  in  the  simulation.  FFT  end 
effects  are  avoided  by  requiring  the  number  of  points  in  a  chip  and  the  number  of  chips  both  to  be  a 
power  of  two  [3].  With  GPS,  the  carrier  frequency  is  an  integer  multiple  of  the  code  frequency,  which 
guarantees  an  integer  number  of  carrier  cycles  in  the  FFT  data. 


RESULTS 

The  simulation  had  different  filters  in  three  areas:  the  receiver,  the  antenna,  and  the  RF  source.  Each  area 
has  three  identical  filters.  Three  different  bandwidths  were  used  for  each  area.  Five  different  values  of 
correlator  spacing  were  used  as  a  fourth  variable  to  cover  both  C/A  and  P/Y.  All  combinations  of  the  four 
variables  produced  a  total  of  135  simulations. 
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Table  1  shows  the  selection  of  variables  used  and  receiver  calibrations  of  the  simulation  runs.  Figure  5 
shows  how  well  different  receivers  track  with  changes  in  simulator  filters.  Figure  6  shows  how 
dependent  receiver  calibration  is  on  antenna  filters. 


CONCLUSIONS 

The  conclusions  are  based  on  the  results  of  the  135  simulations.  Changes  in  the  bandwidth  of  the  RF 
source  will  track  to  within  less  than  2  ns.  Calibration  of  the  receiver  in  parts  can  lead  to  errors  of  6  ns. 
These  are  extreme  condition  and  in  most  cases  they  would  not  be  encountered. 

Since  the  bandwidth  of  the  RF  source  is  typically  greater  than  the  antenna  system,  which  is  typically 
greater  than  the  receiver  bandwidth,  information  on  the  receiver  filters  is  the  most  critical  in  this  analysis. 
Bandwidths  external  to  the  receiver  become  more  significant  as  they  approach  the  receiver  bandwidth. 
Accurate  calibration  of  receivers  with  wide  bandwidths  and  narrow  correlators  is  more  dependent  on 
knowledge  of  the  hardware  than  receivers  with  wide  correlators  and  narrow  bandwidths. 

To  accurately  calibrate  a  receiver  for  time  transfer,  the  hardware  filters  could  be  measured  on  the  network 
analyzer  and  this  information  along  with  correlator  spacing  could  be  used  in  the  simulation.  The  results 
of  the  simulation  then  could  be  used  to  modify  the  receiver  calibration.  An  estimate  of  the  filters  could  be 
used  to  reduce  the  calibration  errors  if  it  is  not  practical  to  measure  the  filters.  With  knowledge  of  the 
hardware,  the  limiting  factor  of  the  receiver  calibration  method  is  the  calibration  of  the  GPS  signal 
simulator  at  the  nanosecond  level  [11- 
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Figure  1.  Block  diagram  of  receiver  calibration  method. 


File  Control  Setup  Measure  Utilities  Help 


Figure  2.  Zero  crossing  of  single  PRN  L  band  simulator  signal  seen  on  a  digital  scope. 
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Figure  3.  Single  pole  filter  envelope  and  relative  position  of  three  measurements. 


Figure  4.  Block  diagram  of  simulation  processing. 
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TIMING  CHANGES  WITH  RF  SOURCE  BANDWIDTH  CHANGE 


—  100  to  50 

—  200  to  50 
200  to  100 


Figure  5.  Changes  in  receiver  calibration  with  changes  in  source  bandwidth. 


TIMING  CHANGES  WITH  ANTENNA  BANDWIDTH  CHANGE 


COMBINATIONS  of  HARDWARE 


—  50  to  none 
-*-10  to  none 
10  to  50 


Figure  6.  Changes  in  receiver  calibration  with  changes  in  antenna  bandwidth. 
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Table  1.  Simulation  bandwidths  and  results. 


Receiver  Calibration  ns 


Antenna 

BW  MHz 

10 

10 

10 

50 

50 

50 

OO 

OO 

OO 

Source  BW 
MHz 

50 

100 

200 

50 

100 

200 

50 

100 

200 

Receiver 

Receiver 

BW  MHz 

Spacing  ns 

2 

977 

465.76 

464.80 

464.38 

466.26 

465.24 

464.83 

466.27 

465.31 

464.89 

2 

313 

435.33 

434.25 

433.84 

433.92 

432.85 

432.43 

433.87 

432.79 

432.32 

2 

97.7 

430.20 

429.12 

428.65 

428.37 

427.30 

426.88 

428.26 

427.19 

426.77 

2 

31.3 

429.66 

428.59 

428.17 

427.84 

426.70 

426.28 

427.72 

426.65 

426.17 

2 

9.7 

429.60 

428.53 

428.11 

427.78 

426.64 

426.23 

427.66 

426.59 

426.1 

10 

977 

97.46 

96.44 

96.02 

97.48 

96.47 

96.05 

97.49 

96.47 

96.06 

10 

313 

95.37 

94.35 

93.94 

96.83 

95.81 

95.39 

96.83 

95.88 

95.46 

10 

97.7 

88.51 

87.37 

86.96 

90.26 

89.13 

88.71 

90.15 

89.02 

88.54 

10 

31.3 

87.31 

86.18 

85.70 

88.12 

86.86 

86.39 

87.76 

86.51 

85.97 

10 

9.7 

87.14 

86.00 

85.58 

87.88 

86.62 

86.09 

87.53 

86.21 

85.74 

20 

977 

49.73 

48.71 

48.29 

49.69 

48.68 

48.26 

49.70 

48.68 

48.27 

20 

313 

48.89 

47.88 

47.46 

49.69 

48.68 

48.26 

49.70 

48.68 

48.27 

20 

97.7 

42.57 

41.49 

41.02 

48.02 

47.13 

46.77 

48.21 

47.31 

46.95 

20 

31.3 

40.96 

39.76 

39.29 

45.75 

44.62 

44.08 

45.52 

44.27 

43.73 

20 

9.7 

40.78 

39.58 

39.11 

45.40 

44.20 

43.67 

45.11 

43.79 

43.12 
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Questions  and  Answers 

DEME TRIOS  MATSAKIS  (U.S.  Naval  Observatory):  Those  are  pretty  large  changes  in  the  receivers 
and  transmitter  systems.  If  you  had  milder  changes  such  as  you  might  get  with  a  temperature  variation, 
calibrating  at  one  temperature,  and  then  installing  at  another,  would  you  expect  it  to  be  a  lot  less 
numerically? 

PAUL  LANDIS:  Yes. 

DAVE  HOWE  (National  Institute  of  Standards  and  Technology):  Can  you  comment  about  what  kind 
of  temperature  stability  you  might  need  for  the  filter,  for  example?  Because  this  also  has  relevance  to 
two-way. 

LANDIS:  I  do  not  have  a  good  feel  for  how  filters  change  with  temperature. 
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